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Band structure calculations have been used to identify the different bands contribut- 
ing to the polarisation-dependent photoemission spectra of the undoped model cuprate 
Sr2Cu02Cl2 at the high-symmetry points of the Cu02 plane F, (tt/o, 0) and (7r/a,7r/a) 
and along the high-symmetry directions T — (n/a,n/a) and V — (n/a,0). Results from 
calculations within the local density approximation (LDA) have been compared with cal- 
culations taking into account the strong electron correlations by LDA+U, with the result 
that the experimental order of energy levels at the high-symmetry points is better de- 
scribed by the LDA+U calculation than by the simple LDA. All the main peaks in the 
photoemission spectra at the high symmetry points could be assigned to different Cu 
3d and O 2p orbitals which we have classified according to their point symmetries. The 
dispersions along the high-symmetry directions were compared with an 11-band tight- 
binding model which was fitted both to the LDA+U band structure calculation and the 
angle-resolved photoemission data. The mean field treatment successfully describes the 
oxygen derived bands but shows discrepancies for the copper ones. 



PACS-numbers: 74.25 Jb, 74.72 Jt, 79.60 Bm 



I. INTRODUCTION 



One strategy to answer the many questions concerning the electronic structure of cuprate superconduc- 
tors is the study of model substances. One of these compounds is r Sr2Cu02Cl2. It is a two dimensional 
(2D) antiferromagnetic (AFM) insulator with a Neel temperatureEl of 256 K whose magnetic structure 
is well described by the 2D spin 1/2 Heisenberg-model. R was the first undoped cuprate which allowed 
the angle resolved photoemission (ARPES) measurement of its lowost excitations.l3a These excitations 
are well described by one hole in a 2D quantum antiferromagnetlftl Deviations from the one- hole dis- 
persion ofrthe pure t-J model can be reduced by taking into account hopping terms to second and third 
neighbors £ru In the same substance the low binding energy edge of the main valence band has been 
interpreted in terms of non-bonding oxygen orbitals which ace completely decoupled from the copper 
system]!!] (These features were known before as "1 eV-peak" These non-bonding states are especially 
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pronounced at (tt, tt) where they have minimal binding energy. But the detailed structure of the complete 
valence band has never been analyzed up to now and that is the aim of the present work. Furthermore, 
we will show that one can obtain additional information on the low binding-energy features by analyzing 
their dependence on the polarization of the photon. 

Polarization dependent photoemission measurements are an effective tool to analyze the electronic 
structure of the valence band in detail. By measuring along high-symmetry directions all bands can 
be classified according to their symmetry properties. This allows a very precise comparison between 
experiment and different theoretical predictions. It is well established that in all the cuprates electron 
correlations have a strong influence on the electronic bands near the Fermi level which is especially 
pronounced in undoped substances B But the influence of correlations on those parts of the valence band 
with larger binding energies is less clear. We will show that the combination of polarization dependent 
ARPES measurements with theoretical investigations taking into account the electron correlation to a 
differing extent (LDA, LDA+U) provides a unique possibility to answer this question. 

The model cuprate S^CuO^CL) is xety well suited for such an investigation. It has a tetragonal 
structure with ideal planar CuC>2 layersEj and cleaves readily parallel to the CuC>2 planes. Furthermore, 
the presence of CI instead of apex oxygen allows a restriction of the states which contribute to the ARPES 
spectra to those of the pure Cu02 plane alone. This can be achieved by choosing a photon energy close 
to the Cooper minimum for CI 3p photoemission, i.e. CI 3p states will then have a small photon cross 
section. In this manner we intend to study an ideal situation whose main characteristics should be generic 
to all the cuprates. |— , 

Recently, a similar study was presented for the CU3O4 plane of Ba2Cu3O.4Cl2.li3 It turns out that 
Sr2Cu02Cl2 is considerably less complex than Ba2Cu304Cl2- Therefore, we are now able to identify all 
the peaks at the high symmetry points in contrast to Ba2Cu304Cl2 where only the upper parts of the 
valence band were analyzed. This allows implications about the influence of electron correlations on the 
valence band (VB) structure of Sr2Cu02Cl2 to be. made. It is known that simple LDA fails to predict 
the insulating ground state of undoped cuprates ££1 These are several improvements of LDA such as the 
self-interaction correction (SIC) methodt-3 or LDA+U£3 which has already been applied to the case of 
lanthanum cuprateU (La2Cu04). Here, we apply LDA+U plus a symmetry analysis at special /c-points 
to interprete the polarization dependent photoemission data for Sr2Cu02Cl2, where the actual value of 
U is chosen to describe the experimental situation. 

The paper is organized as follows. After describing the experimental method and the details of the LDA 
band structure calculation we analyze the symmetry properties of the wave function along high-symmetry 
directions. The symmetry properties of the relevant bands are most clearly seen in a tight-binding model 
presented in Sec. IV. Discussing the correlation effects in a mean-field manner leads us to an LDA+U 
calculation whose results are presented in Sec. V. In Sec. VI we compare the experimental spectra with 
the theoretical predictions. The spectra at T, (71", 0) and (tt, 7t) (the lattice constant has been set to 
unity in all the notations) can be understood from the LDA+U but not from the LDA calculation. The 
experimental dispersion relations are discussed in terms of a tight-binding (TB)-model. 

II. EXPERIMENTAL 

The Sr2Cu02Cl2 single crystals were grown from the melt, their typical dimensions being 3 x 3 x 0.5 
mm. The crystals were mounted on the sample holders using conducting, or in some cases, insulating 
epoxy. If insulating epoxy was used, electrical contact between sample holder and sample was achieved 
by means of a graphite layer at the sides of the crystal. The orientation of the single crystals was 
determined ex-situ by x-ray diffraction. The surface normal of the crystals is perpendicular to the Cu02 
planes. Prior to the ARPES measurements, a clean crystal surface was prepared in ultra high vacuum 
(UHV) by stripping off an adhesive tape which was attached on the sample surface. 

The ARPES measurements were performed using linearly polarized 35 eV photons from the Grossed 
undulator beamline U2 of the BESSY I facility and BESSY's HIRES photoelectron spectrometer M The 
angular resolution was set to ±1° which gives a momentum resolution of ± 0.05 A -1 for states of 1 eV 
binding energy, this corresponds to 12 % of the distance between T and (ir,0). A total energy resolution 
(resulting from both the monochromator and electron analyzer resolutions) of 150 meV was applied. The 
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electron analyzer is fixed in the horizontal plane at an angle of 60° with respect to the incoming photon 
beam, i.e. the emission plane which is defined to be the plane spanned by the sample surface normal and 
the k vectors of the measured photoelectrons, is always a horizontal plane. The photoelectron momentum 
vector could then be changed by variation of the angle between the axis of the input lens of the electron 
analyzer and the sample surface normal by rotation of the manipulator on which the sample holder was 
mounted. The ARPES spectra have been recorded in the electron distribution curve (EDC) mode from 
both r to (tt,tt) and T to (tt,0). In the latter case, the sample was oriented in such a way that the Cu 
- O bonds of the CuC>2 plane were aligned parallel to the fixed, horizontal emission plane, while in the 
former case, the Cu - O bonds were aligned at an angle of 45° with respect to the emission plane. The 
spectra are normalized with respect to the incoming photon flux which was simultaneously measured 
using the drain current of a gold mesh. The absolute binding energy scale was determined according to 
Eb = hv — Ekin — 4>anaiyzer using the photon energy, the measured kinetic energy of the photoelectrons 
and the known analyzer work function. There were no indications of charging-induced energy shifts. All 
spectra were recorded at 300 K within 8 hours of a cleave, during which time the samples showed no 
indications of surface degradation. 

The electronic structure of the VB is derived from O 2p, Cu 3d and CI 3p orbitals, but for 35 eV 
photon energy, the photoionization cross section of the CI 3p orbitals is much smaller than that of the 
O 2p and Cu 3d orbitalsE3, which therefore dominate the ARPES VB spectra. The measurements had 
been performed at room temperature which is slightly above the Neel temperature of Sr2Cu02Cl2 (256 
K). Although we are aware that AFM fluctuations are important, we nevertheless analyzed the spectra 
in terms of the first Brillouin zone (BZ) of the paramagnetic Cu02 plane of Sr2Cu02Cl2. We will see 
that this is especially justified for the bands with dominant oxygen character, whereas one observes 
some deviations for those bands which couple strongly with the copper spins. The geometrical structure 
of a Cu02 plane has two mirror planes (denoted Mi and M2 in Fig. 1)EJ All the bands with a wave 
vector between T and (tt,tt), parallel to the mirror plane Mi, can be classified to be either symmetric or 
antisymmetric with respect to Mi, and analogously also for wave vectors along T — (7r,0) with respect 
to reflections at M2. Experimental information about the parity of the valence band states with respect 
to a mirror plane can be obtained by recording the ARPES spectra with either perpendicular or parallel 
polarization of the electric field vector of the incoming radiation with, respect to an emission plane which 
is parallel to a mirror plane of the system. It can then be shownEil that for parallel polarization only 
initial valence band states which are even with respect to the emission/mirror plane contribute to an 
ARPES spectrum while for perpendicular polarization, only states which are odd with respect to the 
emission/mirror plane are seen in a spectrum. In this work, the polarization of the electric field vector 
was chosen by using either the vertical or horizontal undulator, which corresponds to perpendicular and 
parallel polarization with respect to the emission plane. The emission plane is parallel to mirror plane 
Mi, if the ARPES spectra are recorded along the (7r,7r) direction, while it is parallel to the mirror plane 
M2 for spectra along the (7r,0) direction. For perpendicular polarization, the electric field vector is always 
parallel to the Cu02 planes, i.e. only in-plane orbitals as O 2p x<y or Cu 3d x 2_ y 2 contribute to the spectra. 
For parallel polarization, the electric field vector is completely in-plane only at normal-incidence, at any 
other incidence angle of the photon beam, the electric field vector has an out-of-plane component and 
there are also contributions from out-of-plane orbitals such as O 2p z to the ARPES spectra. 



III. BAND STRUCTURE CALCULATIONS 



Sr2Cu02CL4ias a body centered tetragonal crystal structure with the lattice constants a — 3.973 A and 
c = 15.618 A.li3 Band structure calculations have been performed treating the exchange and correlation 
potential within the local density approximation (LDA). The Bloch wave functions were constructed from 
atomic-like wavefunctions according to the linear combination of atomic orbitals (LCAO) method. The 
calculation was scalar relativistic and due to the open crystal structure two empty spheres per elementary 
cell were introduced in between two oxygen atoms of neighboring Cu02 planes. A minimal basis was 
chosen consisting of Sr(5s,5p,4d), Cu(4s,4p,3o?), 0(2s,2p), Cl(3s,3p) orbitals and the Is and 2p orbitals 
for the empty spheres. To optimize the local basis a contraction potential (r/ro) 4 was introduced.EHl The 
Coulomb potential is constructed as a sum of overlapping contributions of spherical symmetry and for 
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the exchange and correlation potential the atomic sphere approximation (ASA) is used. 

In the resulting band structure (Fig. 2) one observes an antibonding band built up of Cu 3d x 2_ y 2 
and O 2p xy orbitals crossing the Fermi level. This contradicts the experimentally observed non-metallic 
behavior which already indicates that one has to treat the electron correlations in a more explicit way. 
One could conjecture that the only effect of correlations is to split the half-filled antibonding band leaving 
the structure of the other valence bands roughly unchanged. That is not the case, however, as will become 
clear from our following analysis. One can also observe in Fig. 2 that there is nearly no dispersion of the 
relevant band in the z direction and all discussions in the present paper will be restricted to the Cu02 
plane only. 

To check the minimal basis LCAO method, a full potential linearized augmented plane wave (FLAPW) 
calculation has been performed for Sr2Cu02Cl2 (Fig. 3, compare also ReL-23). This method involves 
no shape approximations and uses a flexible basis in all regions of space.ESAs such it is well suited 
to open structures with low site symmetries as in the present cuprate. We note a sufficiently good 
agreement between both band structures, although the LCAO-bandwidth of the valence band is found 
to be somewhat larger. However, there are no significant differences in the order of energy levels between 
LCAO-LDA and FLAPW-LDA. We stick to the LCAO because we want to exploit the minimum basis 
orbital analysis. 

To obtain more information about the structure of the valence band in our LCAO-LDA we have 
calculated the orbital weight (defined in Ref. ^5|) of each band at the high symmetry points. Due to the 
low cross section of the CI 3p orbitals for 35 eV photon energy we concentrate on the Cu 3c? and O 2p 
orbitals (i.e. on 11 bands). The eigenfunctions with a dominant contribution from Cu 3d and O 2p orbitals 
are collected in Table 1. The in-plane oxygen orbitals are divided into p a orbitals which are directed to 



the Cu site and p w orbitals perpendicular to themES There are two combinations for each: p a and p a , 
(Pk and Pk), which are antisymmetric and symmetric with respect to reflection in Mi, respectively. The 
precise definition of these orbitals will be given in the next section. 

Thus we are able to predict the symmetry of each band at the high symmetry points in the Brillouin 
zone (BZ). However, as will be seen later, the order of energy levels of the LDA calculation is incompatible 
with the experimental spectra. Moreover, as it was mentioned already, LDA calculations are unable to 
describe the Mott insulating ground states of the undoped cuprates and do not produce the Cu local 
moments that are present in these systems. The splitting of the spectral density due to the 3d x 2_ y 2 
states away from the Fermi energy due to Coulomb correlations and the resulting reduction in Cu-0 
hybridization is expected to be largely missing in such calculations. However, what, if any, changes 
there are from the LDA bands away from Ep is unclear, particularly well above T/v, where the magnetic 
scattering due to antifcrromagnetic spin fluctuations should be more or less incoherent. Addressing this 
question is one of the main goals of the present paper. In the following we develop a more sophisticated 
LDA+U calculation taking into account explicitly the effects of strong correlations. As a preliminary step 
we formulate an effective tight-binding model which will be fitted both to the LDA+U band structure 
calculations and the ARPES VB spectra. 



IV. SYMMETRY ANALYSIS AND TIGHT-BINDING MODEL 



The polarization dependent ARPES measurements of VB states along the two high-symmetry direc- 
tions r — (jr, 7r) and r — (-7T, 0) discriminate the parity of these states with respect to reflections in the 
corresponding mirror planes Mi and M2. To make the analysis of the experimental data more straight- 
forward it is helpful to incorporate the symmetry properties of the VB states in our approach from the 
beginning. This becomes especially clear by constructing an effective tight-binding (TB) model taking 
into account the point-group symmetry of the VB states. The TB model will be restricted to the 11 
bands of Cu 3d and O 2p. Of course, as can be seen in Table 1, there occurs in some cases quite a strong 
mixing with the CI subsystem, but in the following we will assume that this mixing is taken into account 
by the particular values of the TB parameters. 

We start with the description of in-plane oxygen orbitals whose analysis is more involved than that 
for the copper or out-of-plane orbitals. We introduce the annihilation operator of an electron in the two 
oxygen 7r-orbitals belonging to an elementary cell at position i (i is a site of the square lattice) as P^g/ 2 t 
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where (a, (3) — (x, y) or (y, x) with x and y to be the two orthogonal unit vectors of the lattice. The 
d xy orbitals hybridize with a particular combination of oxygen orbitals arranged over the plaquette at 

site i: p^i = 2"(Pi+L/2 ~ Pi- x /2 ~^Pi+ y /2 ~ P^-y/2)- T ms plaquette's 7r-orbitals are not orthogonal to each 
other. The orthogonalization can be made by introducing first the Fourier transformation for the original 
p^-orbitals 

~ JN ^ 1+0,2 

i 

At the second step we define two kinds of canonical Fermi-operators 

P „{ q ) = \-\{ Sq j*\ q )- Sq ^:\ q )) 

Ml) = K^'P^b) + s g ,yPi x \q)) (1) 

where s 9iQ = sin(q Q ,/2) (a — x,y) and A g = 1 s 2 , x + s 2 , ~. It is easy to see that p n and p„ are orthogonal 
with respect to each other. The definition (1) provides an equivalent representation for 7r-orbitals in terms 

of p v (q) and p n (q) , instead of the original (<?) and p^f 1 (q) operators and takes into account the point 
group symmetry of the Cu02 plane. In particular, for q along T — (w, it), the p^-orbital is antisymmetric 
with respect to reflections in the mirror plane Mi, while the p^-orbital is symmetric (see Fig. 4). Along 
r — (tt, 0), we find p n to be symmetric and p n to be antisymmetric with respect to reflection in M2. 

Turning now to the oxygen er-orbitals we carry out the same procedure as above with the corresponding 
Pi+a/2 a °P era tors (a = x,y). In this case, introducing the plaquette representation instead of defining 
the original p„ operators in momentum space, we define a new pair of canonical Fermi-operators p a and 

Pa- 



P<t(q) = \ li ( s q,yPq X J + Sq,xP q V J) ■ 

The notation is chosen in such a way that the p CT ((7)(p CT (<7))-orbitals have the same symmetry properties 
with respect to reflections at Mi and M 2 as the p n (q) or p T (g)-orbitals, respectively. 

The definition of the corresponding copper annihilation operators is quite standard and thus we may 
write down the TB Hamiltonian 

H t =Y,cUl)H^{q)c vs {q). (2) 

q^ivs 

Here, c MS is an annihilation operator of either an oxygen p orbital or a copper d orbital, where the indices 
H and v denote the 11 different orbitals and s denotes the spin. All orbitals can be classified as to 
whether they hybridize in-plane or out-of-plane and there is no coupling between the two subsystems. 
The orbitals involved in the hybridisation in-plane are p a , p v , p a , p w , d x 2_ y 2, d xy , d^ z 2_ r 2. The explicit 
form of the TB-Hamiltonian for in-plane orbitals is given in the Appendix. The in-plane part of the TB- 
model has 11 parameters: the on-site energies Ed (for d x 2_ y 2), Ed (for d xy ) and (for d 3z 2_ r 2) as well as 
e p (corresponding to p a ) and e^; the hopping matrix elements t p d, t p D, t j, t pp , and t p7T . Besides the 
orbitals hybridising in- plane we have to consider those involved in hybridization out-of-plane: O 2p z , Cu 
3d xz and Cu 3d yz . Restricting ourselves to nearest neighbor hopping leads to two 2x2 matrices with 
on-site energies £ pz and Edz and the hopping matrix element t p d z . 

In order to analyze the experiment it is important to know the parity of the orbitals with respect 
to reflections at the corresponding mirror planes Mi and M 2 . This can also be expressed in terms of 
group theory since for k vectors along the line Y — Ir^^f 1 ) au wave functions can be classified in terms 
of irreducible representations of the small group C^v E3cj The bands built up from the in-plane orbitals 
d xy , d 3z 2_ r 2, p a and p^ belong to the representation Ai and are symmetric with respect to reflections at 
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Mi, whereas d x 2_ y 2 : p a and p v belong to A2 and are antisymmetric. The same small group Civ also acts 
along T— (71", 0) and the subdivision of the in- plane orbitals is as follows: A\ (symmetric): d x z_ v 2 , d 3z 2_ r 2 , 
p a , and A2 (antisymmetric): d xy , p„, p^. The subdivision along high-symmetry lines is also easily seen 
in the TB-matrix given in the Appendix. The corresponding small groups at the high-symmetry points 
r, (?r, 7r) and (tt, 0) are and D 2 h, respectively, and the assignment of the different orbitals to the 
corresponding irreducible representations is given in Table 2. Of course, the group theoretical analysis is 
not only valid for the TB model but also for the LDA bands (Table 1). 
The TB-Hamiltonian (||) should be completed by an interaction term 

H = H t + Hu (3) 

which will not be written out explicitly. This is just a direct extension of the three-band Emery model to 
the case of the complete set of 11 bands for the Cu02 plane. The interaction term Hu involves intrasite 
Hubbard repulsion for different kinds of copper and oxygen orbitals and appropriate intersite copper- 
oxygen repulsions. In order to establish the one-electron parameters entering into Eq. ([}]) one has to 
keep in mind that these parameters are "bare" ones while the results of the band structure calculations 
should be interpreted in terms of a mean-field solution of Eq. (||).l!j To arrive at the bare parameters, one 
would have to take into account the ground-state (GS) properties of the Cu02 plane and approximate 
the Coulomb interaction terms. 

In an undoped cuprate compound as S^CuO^C^, the GS of a particular Cu02 plane contains one 
hole per cell which is shared between d x 2_ y 2 and p a orbitals. Thus a convenient description of the GS is 
to introduce the deviations (n d )h = 1 — (n d ) and {n p )h = 1 — (n§) from the full band (Cu 3d 10 O 2p 6 ) 
electron occupancy. A rough estimate is (n d )h ~ 0.7 and (n p )h ~ 0.3. Here (n p ) means the electron 
number in the p a orbital with spin s (one should note that the occupation of a local oxygen orbital is 
only half that number). Now the mean-field ("screened") one-electron energies £ MS read as follows 

e ds = e d + U d - U d (n d s ) h - 2U pd (^2 n p s ')h 

s' 

£ ps =e p + U p - 2U pd (^2 n t>)h - ^U p (n p g ) h 

s' 

e D = e D + U d - U dD (J2n d s ')h - 2U Dp {J2n p ,) h 

s' s' 

e 7r = £ p + U„- -U^C^n^u - 21/^(^2 n d s ,) h , (4) 

s' s' 

where s — —s. There are also similar expressions for ej, E dz , e pz which we do not specify here. 

In the paramagnetic LDA band structure where the correlation effects are treated only in an averaged 
manner, the screening effect is nearly the same for all ci-levels. So, in the LDA approach the effects 
of strong correlations due to U d are missed. An obvious way to adopt these effects is to treat the 
ferromagnetic solution by putting, for instance, (n d )h = 0, and {ni)h = n d . Then e d ] — e d + U d — 2U pd n p 
(nP = J2 s ( n ^) h )> ^ s lifted upwards while e dy = e d + U d (l — n d ) — 2U pd n p is shifted equally downwards 
with respect to the paramagnetic solution. Regarding the other d-levels, let us assume for the moment 
the rough estimate for the intrasite Coulomb parameters U d o — U d . Then one can see that £d = 
£d + U d (l — n d ) — 2U pd n p , and the d xy level as well as all the other remaining Cu (i-levels are shifted 
as was the lower e^. The spin dependence of e ps in (Q) is much less pronounced than for e ds and is 
neglected in the following. 

Thus, although being somewhat awkward, the ferromagnetic solution provides a better description of 
the strong electron correlations, giving a more reasonable energy position and occupancy of the different 
orbitals. Just this approach is taken by us to carry out the LDA+U calculation. The details of the 
procedure and some results of these calculation are presented in the next Section. 
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V. LDA+U CALCULATION 



The main effect of a mean-field treatment of the multi-band Hubbard model is a shift of the on-site 
copper energies against the oxygen ones. Furthermore, the on-site energy of the Cu 3d x 2_ y 2 orbital is 
split into one for spia up (minority spin) and Sd[ (majority spin). This can also be achieved by an 
LDA+U calculation^ including all valence orbitals. 

We performed LDA+U calculations for Sr 2 Cu02Cl2 using a ferromagnetic splitting. The on-site energy 
of the unoccupied, spin up Cu 3d x2 _ y2 orbital (minority spin) is shifted by 2 eV upwards and the occupied, 
spin down Cu 3d^ 2 _ i) 2 orbital (majority spin), as well as both spin directions for all the remaining Cu 3d 
orbitals are shifted by 2 eV downwards. The energy shifts were added at each step of the self-consistency 
cycle until the charge-distribution was stable. We did not try to connect the chosen energy shifts with the 
model parameters such as, for instance, Ud,Upd,U p . According to (|4|), the actual shift depends also on 
the occupation numbers {nf)h and {n p s )h- Since we did not shift the oxygen levels, our choice corresponds 
in fact to the difference between Ud and U p weighted with the corresponding occupation numbers. 

The results of our LDA+U calculation are presented in Fig. 5 and Table 3. The mainly unoccupied, 
minority band of d x 2_ y 2 and p a character can be roughly interpreted as the upper Hubbard band. The 
corresponding band for majority spin lies just below the Fermi level and has dominantly oxygen character. 
Since its spin is opposite to the spin of the copper hole, there is some justification to interprete that band 
as the mean field representation of the Zhang-_B,ice singlet. But due to our ferromagnetic spin structure it 
has completely the wrong dispersion relation.lEl The bandwidth of both bands is expected to be strongly 
reduced by correlation effects in comparison with Fig. 5 such that a gap opens. 

Next in binding energy we find bands with dominantly oxygen character. The nonbonding oxygen band 
with lowest binding energy at (n, ir) is identified to be of pure p^ character. The oxygen bands occur at 
nearly the same energy for both spin directions. In fact, only the bands with a considerable weight of 
the Cu 3d x 2_ y 2 orbital show a strong splitting between spin up and spin down. Therefore we present in 
Table 3 only the position minority spin bands and both spin directions for bands with a contribution 
from the 3d x 2_ y 2 orbital.tj The actual value of the energy shifts of the copper bands in our LDA+U 
calculation has little influence on the upper oxygen bands, only their copper character is changed. We 
have chosen such a shift that the copper bands are at the lower edge of the valence band, but are not yet 
split off the valence band. This is important in order to achieve good agreement with the experimental 
results. 

Let us now compare the LDA and LDA+U results starting at (7r,7r). In both cases (Figs. 2 and 5), 
we find a group of 5 bands at around 3 eV binding energy, but the order of energy levels is completely 
different in the two cases. For example, the antisymmetric p^ band has lowest binding energy of f=a 2.5 
eV in the LDA+U calculation. In Fig. 2 (LCAO-LDA), however, all the other 4 bands of that group 
have lower binding energy than the p„ level. And also in the FLAPW calculation (Fig. 3) the p^ band 
has 0.3 eV larger binding energy than the valence band edge. A similar rearrangement of energy levels 
can be observed at the V point. Due to symmetry reasons there is no hybridization between copper and 
oxygen bands there. The energy position of the oxygen bands is nearly the same for LDA and LDA+U, 
but the copper bands are shifted. The in-plane oxygen bands are twofold degenerate and occur twice in 
the LDA+U result with binding energies of 2.69 and 5.57 eV, respectively. 

VI. COMPARISON WITH EXPERIMENT 

A. High-symmetry points 

The experimental ARPES spectra at the high symmetry points for both polarization directions are 
presented in Fig. 6. At the T point, there are two possible orientations of the sample such that one 
can probe the symmetry of states with respect to reflections in either Mi (sample directed such that the 
photoelectron momentum is along T—(ir, ir), Fig. 6a), or M2 (sample directed such that the photoelectron 
momentum is along T — (%, 0), Fig. 6b). The first peak at 2.9 eV binding energy in the experimental 
spectra at the T point with the sample oriented such that the fc-vector is along T — (ir, tt) (Fig. 6a) is 
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equally strong for both polarization directions. This leads us to interprete it as the two pure oxygen 
bands (p^paX, and (p v pa) which are antisymmetric and symmetric with respect to reflections at Mi, 
respectively.t3 These bands occur in the LDA+U calculation as the two-fold degenerate in-plane oxygen 
bands at 2.69 eV binding energy. According to this interpretation we would expect the same identical 
peak for both spin directions also at the T point with the sample oriented such that the fc-vector is 
along r — (tt, 0) (Fig. 6b). As one can see, Fig. 6b deviates only slightly from that expectation. In the 
LDA result, however, there are 3 copper levels between 2.3 and 3 eV binding energy. Since every copper 
level has different symmetry properties with respect to Mi and M2 that would lead to strong differences 
between both polarization directions which is not observed. Therefore, we assign each experimental peak 
with the help of the LDA+U results. Each pure band is denoted by one orbital only. For the mixed bands 
we choose a notation using two orbitals, where the first one is the dominant one. The experimental peak 
positions are compared with the LDA+U positions in Table 4. 

Let us continue our interpretation of the spectra at the T point with the peak at 3.9 eV. It is seen 
with horizontal polarization in Figs. 6a and 6b. Therefore, we interprete it as the out-of-plane oxygen p z 
orbital. We observe also a small contribution of this peak with the "wrong" polarization in Fig. 6a which 
is even larger in Fig. 6b. However, there is no band with the corresponding symmetry in that energy 
region in our LDA+U calculation. The large peaks at around 6 eV binding energy in Figs. 6a and 6b 
with big differences between both polarization directions indicate that there are additional contributions 
besides the oxygen orbitals there. Due to the low cross section of CI 3p orbitals, we are only left with the 
pure copper d orbitals. To simplify the analysis we did not try to assign the Cu 3d 3z 2_ r 2 orbital which 
mixes strongly with the CI orbitals and should have reduced intensity. The remaining in-plane copper 
orbitals change their polarization dependence between Fig. 6a and 6b. The d 2 _ 2 is antisymmetric with 

x y 

respect to Mi and the d xy is symmetric, but with the sample oriented such that the fc-vector is along 
L — (tt, 0) this situation is reversed. The intensity ratio between horizontal and vertical polarization of 
the peak at 5.8 eV is indeed exchanged if we compare Fig. 6a and 6b. The last peak at 6.5 eV occurs for 
both sample orientations only with horizontal polarization and is interpreted as the out-of-plane d xz or 
d yz orbital. 

Turning now to the spectra at (tt,tt) we can only probe the parity with respect to Mi (Fig. 6c). The 
small piijepeak at 1.2 eV in the curve with vertical polarization is ususally interpreted as the Zhang- Rice 
singlet.B The dominant peak at 2.4 eV binding energy in the spectra with perpendicular polarization 
can be identified as the pure p^ orbital which has already been discussed in Ref. [l0|. The p^ band is 
the only one among the group of 5 bands at around 3 eV binding energy in both calculations (LDA or 
LDA+U, Figs. 2 and 5) which has odd symmetry with respect to Mi. It has lowest binding energy in 
the experiment and in the LDA+U calculation. That indicates that the LDA+U calculation is better in 
predicting the correct order of energy levels at high symmetry points than the pure LDA calculation. At 
slightly higher binding energy at 2.7 eV we observe a smaller, broader peak with horizontal polarization. 
According to our calculation it should be comprised of three bands, the out-of-plane (p z d( X4/ - )z ) bands 
and the in-plane (p-„d X y) band. The small structure at 3.8 eV binding energy (vertical polarization) can 
be related to the oxygen p a orbital hybridizing with d x 2_ y 2 but having opposite spin (|) than that of 
the copper hole. The corresponding band occurs in the LDA+U at 4.94 eV binding energy and can be 
interpreted as the Zhang-Rice triplet. A similar structure was also observed in our previous analysis-pf 
the polarization dependent photoemission spectra of another undoped model cuprate Ba2Cu304Cl2Ji3 

The peaks at around 6 eV binding energy should be assigned to bands with a dominant copper character. 
But one may note in Table 4 a systematic deviation between experimental and theoretical peak positions 
at (tt,tt): the theoretical binding energies are too large. That is plausible since it is expected that 
the copper bands feel the antiferromagnetic correlations much more than the oxygen bands which are 
decoupled from the copper spins. As a result the copper bands are expected to follow more the AFM BZ 
where T and (tt, tt) are identical. However, such AFM correlations were not considered in our calculation. 

At (tt, 0) (Fig. 6d) one may observe a prepeak with low intensity which may be prescribed to the 
Zhang-Rice singlet state comprised in our calculation by the hybridization between the p a orbital and 
d x 2_ y 2- The strong peak with horizontal polarization at 2.5 eV is assigned to the out-of-plane (p z d xz ) 
orbital. The peak at 3.8 eV consists of two orbitals p z and p n which are separated by only 0.5 eV in 
the LDA+U calculation. Therefore it is difficult to use that peak to extract the parameter from the 
experimental spectra as it was done in Ref. fill. Furthermore, one should distinguish between different 
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oxygen hopping matrix elements (t ppt t p7T and jtli which was also not done there.H9 

B. Dispersion relations 

The experimental spectra along both high symmetry directions show clear differences between both 
polarization directions (Figs. 7-10). The first electron removal peak along T — (jr, tt) has minimal binding 
energy at (7r/2, 7t/2) and occurs exclusively with vertical polarization (Fig. 7). That is in complete agree- 
ment with the usual interpretation of that peak as the Zhang- Rice singlet. In our mean-field treatment 
it is built up of the d^ 2 2 and p a orbitals having odd symmetry with respect to Mi. The dispersion 

is well described within the extended t-J modeltl and we have included the corresponding theoretical 
curve in Fig. 7 for completeness. Along T — (tt,0) (Figs. 9 and 10) the Zhang- Rice singlet feature is less 
pronounced and according to our symmetry analysis based on a simple mean-field treatment we would 
expect it only with horizontal polarization. However it is more clearly seen in Fig. 9 (vertical polarisation) 
than in Fig. 10 (horizontal polarisation). The explanation of that effect deserves obviously a more refined 
treatment and will be studied both theoretically and experimentally in the future. 

The peak next in binding energy in Fig. 7 was already analyzed as the p^ orbital and it has a clear 
dispersion going from T to (tt, tt). The valence band edge at around 2.5 eV binding energy is different 
for both polarizations along T — (tt, 0) as well: it has no dispersion for vertical polarization (Fig. 9) and 
is built up of only one (p^) orbital. In contrast to that, we see for horizontal polarization (Fig. 10) one 
dispersionless out-of-plane band at 3.9 eV and two crossing bands from the out-of-plane orbitals and the 
in-plane p^ band. 

To analyse this dispersion quantitatively it is more convenient to use the TB model than the LDA+U 
calculation due to the restricted number of bands in the former. The parameters of the TB model 
were found as follows. The LDA+U results at high symmetry points (Table 3) were used to obtain 
a first parameter set. For the fit we have only chosen such energy levels which have no or very small 
contribution from other orbitals (Cu 4s, O 3s, CI). In such a way our effective TB parameters also contain 
the influence of hybridization to CI or s orbitals. Fitting to the pure LDA results (Table 1) gave nearly 
the same hopping integrals but different on-site energies. The parameters are very similar to those known 
for La2Cu04O After fitting to the LDA+U results there remained small differences to the experimental 
dispersions even for the peaks with lowest binding energy. These small discrepancies to the experimental 
peak positions were corrected by small changes of the on-site and off-site energies (here, especially t p d z 
was increased). The resulting parameter set is shown in Table 5. 

In Fig. 11 we have collected all the peak positions from Figs. 7-10 together with the dispersion of the 
TB bands. We have distinguished between the results for vertical polarization (Fig. 11a) and horizontal 
polarization (Fig. lib). According to our previous analysis, the peaks in Fig. 11a between (7T,0) and T 
should only be compared with the 3 TB bands stemming from the p„, p^ and d xy orbitals. Analogously, 
between V and (tt,tt) (Fig. 11a) we present only the antisymmetric bands from the p a , p^, d^. 2 _ y2 and 
d} 2 2 orbitals. In Fig. 11, we have collected the bands arising from both the d^ 2 2 or d^ 2 2 orbitals, 

x — y ° ' ° x — y x — y ' 

and have neglected the band corresponding to the Zhang-Rice singlet since we cannot expect to obtain 
its correct dispersion in our simple mean-field treatment. The number of bands which contribute to the 
spectra for horizontal polarization (Fig. lib) is considerably larger: these include all of the out-of-plane 
orbitals and additionally the corresponding symmetric bands (representation A\ of C2 V ) of the in-plane 
orbitals. 

In Fig. 1 1 we have distinguished between bands with dominant oxygen character for all k values (solid 
lines) and those bands which have a considerable coupling to the copper spins (dashed lines). As one may 
note, there is a considerable agreement between experimental and theoretical dispersions for the oxygen 
bands with small binding energy. Furthermore, there is some similarity at the T point besides the peak 
with vertical polarization at 3.9 eV binding energy for which we have no explanation. But the copper 
bands at around 6 eV disperse strongly in the TB calculation whereas they are nearly dispersionless 
in the experiment. We think that this failure of the theoretical description arises due to the neglect of 
antiferromagnetic correlations. To avoid misunderstanding we should stress that also the oxygen bands of 
our mean-field calculation have a copper contribution (except some cases at the high-symmetry points), 
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but that the copper contribution is not dominant. We have also shown the calculated dispersion relations 
of the oxygen bands in Figs. 7 - 10 as solid lines in order to guide the eye. 

VII. CONCLUSIONS 

It can be summarized that polarization dependent ARPES at T, (it, it) and (it, 0) and along the two 
high-symmetry directions gives detailed information about the bands with different parity with respect 
to reflections at the mirror planes Mi and M2. The assignment of the peaks can be performed by means 
of a symmetry analysis of band structure results. Here we pick out the three major results. 

Rearrangement of energy levels. Comparing LDA with LDA+U results at high-symmetry points 
we found that the strong electron correlation leads to a changed order of energy levels, whereby the 
experimental peak positions could be more accurately assigned with the help of the LDA+U calculation. 
In comparison with an LDA calculation we found the copper bands shifted to higher binding energy. 
So, we conclude that the correlation influences not only the band near the Fermi level but leads to a 
rearrangement of energy levels throughout the whole VB. 

Check of the non-bonding p v band. Polarization dependent ARPES measurements provide a 
sensitive test of the symmetries of the excitations with low binding energy which were already analyzed 
before. The p n orbital is seen at (it, it) with vertical polarization as a single peak. At (tt, 0) it is visible 
with horizontal polarization but overlaps with out-of-plane orbitals which makes a parameter assignment 
difficult. This means that in polarization independent measurement, such as those in Ref. the spectral 
weight assigned to the p n peak at (tt, tt) will have additional contributions besides the pure p n orbital (of 
roughly one third of the total intensity as seen in Fig. 6c). As a consequence, the experimental estimate 
of the spectral weight of the Zhang- Rice singlet part, which was performed there using the intensity of 
the Pn feature as a calibration, should be increased by 50 per cent. 

Dispersion relations. Analyzing the dispersion relations we observe a difference beween the copper 
bands which couple strongly to the antiferromagnetic spin structure and thus feel the antiferromagnetic 
BZ and the nonbonding oxygen bands which are decoupled from the spin system and follow the para- 
magnetic (or ferromagnetic) BZ. To take that effect into account for S^CuC^C^ we should extend our 
theory twice. First we should incorporate the antiferromagnetic order. Then all the bands are defined 
within the AFM BZ. To obtain in such a scheme the observed difference between T and (it, it) deserves 
the calculation of matrix elements. 

Despite the fact that the experimental order of energy levels can be explained by an LDA+U calculation 
one should be aware that the agreement between photoemission and LDA+U cannot be perfect. First of 
all, the LDA+U calculation cannot reproduce the satellite structure present in the spectra at about 14 eV 
binding energy. And second, the LDA+U has the tendency to push the copper levels to too large binding 
energy. That was visible in our analysis especially at (it, it). The fc-integrated copper density of states 
can also be measured by x-pty photoemission with large photon energy such that the copper cross section 
dominates that of oxygen.tS It was found that the x-ray photoemission spectrum of the valence band 
of Si'2Cu02Cl2 showed the existence of Cu 3d electron removal states aszer an energy range of some 5-6 
cV. To compare our LDA+U calculation with earlier ones for La2Cu04E3, one should also keep in mind 
that we had to choose a rather small shift of the copper levels to find agreement with the experimental 
situation, and we did not choose correlation parameters from a constrained density functional calculation 
as m Ref. 0. It can be expected that the consideratiaippf self-energy corrections as was done recently 
by calculating the three-body scattering contributionscj'tH improves the situation and allows one to work 
with real correlation parameters instead of fitted ones. Our main goal here was the assignment of peaks 
and not the determination of parameters. To extract parameters from polarization dependent ARPES 
measurements there are several improvements necessary both from the experimental and the theoretical 
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side of view. 
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Appendix 

The TB-matrix for the in-planc orbitals can be found in the form 
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where s q ^ x , s q ^ y and X q arc defined in the main text and the other expressions are given by: 
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FIGURES 



FIG. 1. The mirror planes of the Cu02 plane. Filled (open) circles correspond to copper (oxygen) 
atoms. 

FIG. 2. The LDA-LCAO band structure of Si^CuG^C^. The wave vector is given in units of 
(ir/a, n/a, 7r/c). 

FIG. 3. The LDA-FLAPW band structure of Sr 2 Cu02Cl2. The points in fc-space are denoted as 
Z=(0,0,tt/c) and X=(7r/a,7r/a,0). 

FIG. 4. Scetch of the different oxygen orbitals within one unit cell (filled circles - copper, open circles 
- oxygen) for momenta q — > along T — (tt, tt). 

FIG. 5. LDA+U band structure: (a) minority spin (f), (b) majority spin (J.). 

FIG. 6. Experimental photoemission data at high-symmetry points, with the relevant mirror plane 
given in brackets: (a) at the T point (Mi), (b) at the F point (M 2 ), (c) at (tt, tt) (Mi) and (d) (tt, 0) (M 2 ). 
The assignment of peaks is according to the LDA+U results. The filled circles and full lines correspond to 
vertical polarization, whereas the open circles and broken lines give the results for horizontal polarization. 

FIG. 7. Angle resolved photoemission curves along T — (tt,it) for vertical polarization. Also shown is 
the antisymmetric TB-band with dominant oxygen (p^) contribution (full line) and the dispersion of the 
Zhang- Rice singlet according to Ref. || (dotted line). 

FIG. 8. Angle resolved photoemission curves along T — (tt, tt) for horizontal polarization together with 
the calculated oxygen out-of-plane TB-bands and the in-plane bands having even symmetry (calculated 
bands are shown as solid lines). 

FIG. 9. Angle resolved photoemission curves along T— (tt, 0) for vertical polarization with the calculated 
antisymmetric oxygen TB-band (shown as a solid line). 

FIG. 10. Angle resolved photoemission curves along T — (tt, 0) for horizontal polarization with the 
calculated oxygen out-of-plane TB-bands and the in-plane bands having even symmetry (calculated bands 
are shown as solid lines). 

FIG. 11. Position of the main experimental peaks together with the TB-bands of the correspond- 
ing symmetry along (tt,0) — (0,0) — (tt,tt): (a) antisymmetric bands and experimental data for vertical 
polarization, (b) out-of-plane and symmetric bands together with experimental data for horizontal po- 
larization. Full lines denote the TB bands with dominantly oxygen character, whereas the dashed lines 
correspond to bands with a considerable mixing to the copper system. 
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TABLES 



TABLE 1: LDA data at high symmetry points showing the weights of the different orbital groups 
contributing to each band. Also given are the different reflection symmetries with respect to Mi and M2, 
respectively (antisymmetric (A), symmetric (S) and out-of-plane bands (o)). 
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TABLE 2: Assignment of the orbitals to irreducible representations of the corresponding small groups 
at high symmetry points: a) T (group rQAh ) , b) (ir, 7r) (Djh) and c) (%, 0) {D^hj- The notations in 
parantheses are according to LuehrmannlHa (see also Ref . E6h . Also given are the characters with respect 
to reflections at Mi or M2, respectively, whereby + and - correspond to the S and A given in Table 1. 
The orbital means p z orbitals at positions i ± x/2, and p[ 2 ^ at positions i ± y/2. 



(a) r 



orbitals 


repr. 


Mi 


M 2 


Pa, Pa 


(5~) 








Pn, Pit 


El?> (5-) 










A 2u (2~) 


+ 


+ 




B 2u (4-) 




+ 


dx 2 — y 2 


Big (3 + ) 




+ 




B 2g (4 + ) 


+ 




d(x,y)z 


E g (5+) 








d 3z 2_ r 2 


M 9 (i + ) 


+ 


+ 



(b) (7T,7T) 



orbitals 


repr. 


Mi 


dzz 2 -r 2 , Pa 


M 9 (1 + ) 


+ 


Pit 


A 2g (2+) 




d x 2 —y 2 , Pa 


Big (3+) 




dxy , Pit 


B2g (4+) 


+ 


1 (1,2) 
">(x,y)z, Pz 


Eg (5+) 






(c) (tt,0) 



orbitals 


repr. 


Mi 


d x 2_ y 2, d 3z 2_ r 2, p a 


Ag (1 + ) 


+ 


dxy, Pit 


Big (2+) 




Pa 


B 2u (3-) 




Ptx 


B 3u (4-) 


+ 


Uxz, Pz 


B 2g (3+) 


+ 


(2) 
Pz 




+ 


dyz 


B 3 g (4+) 


+ 



l(i 



TABLE 3: The LDA+U data at the high symmetry points. The bands noted by a star correspond to 
majority spin (J.), whereas all the other data are given for minority spin (j). The column "not." gives 
the notation used to describe the bands. 



r 



No. 


E/eV 


Pz 


Pa 


Ptt 


d 3 z2„ r : 










d x 2_ y 2 


o s 


Cu s 


Eci 


not. 


M a 


M 2 


1 


-.30 























.909 


.091 








d T 2 2 

X Z_yZ 


A 


S 


2,3 


-2.69 





.439 


.543 























.018 


(p^p<t) 

(PttP^) 


S 
A 


A 
S 


4,5 


-3.83 
(-3.69) 


1.00 
(.567) 































(.433) 


Pz 


o 


o 


6 


-4.58 











.640 














.031 





.329 


d 3z 2_ r 2 


S 


S 




-4.92 























.893 


.107 










A 


S 


8 


-5.40 




















1.00 
















S 


A 


9,10 


-5.57 





.499 


.463 























.038 


(PctPtt) 
(PttPo-) 


S 
A 


A 
S 


11,12 


-5.88 














.803 


(.803) 














.197 


<iy(x)z 





o 


(tT, 71") 


No. 


E/cV 


Pz 


P<T 


Ptt 


d 3z 2_ r 2 




^•x z 


da; j, 


d x 2_ y 2 


o s 


Cu s 


Eci 


not. 


M 2 


1 


3.12 





.467 

















.533 











{d l x 2_ y 2Pa) 


A 


2* 


.65 





.704 

















.296 











(P*<&- 


y 2) 


A 


3 


-2.43 








1.00 


























Ptt 


A 


4,5 


-2.98 
(-2.97) 


.711 











.143 


.143 














.002 
(.006) 


(Pzd x (2,)z) 





6 


-3.35 





.009 


.345 


.030 








.564 








.004 


.048 


(PTrd^) 


S 


7 


-3.66 





.096 


.020 


.356 








.053 








.051 


.424 


(d 3z 2_ r 2 


P<x) 


S 


8 


-4.94 





.561 

















.439 











(Pad l x2 _ 


y2) 


A 


9,10 


-6.62 
(-6.06) 


.155 
(.093) 











.341 
(.275) 


.341 
(.275) 














.123 
(.357) 


(dx(y)zPz) 





11 


-7.20 





.003 


.392 


.010 








.591 











.004 


(dxyP7r) 


S 


12* 


-7.28 





.207 

















.793 











d a .2_ w 2Pcr 


A 


13 


-7.86 





.478 


.001 


.073 

















.280 


.168 


(pad 3z 2_ 


-r 2 ) 


S 



(tt,0) 



No. 


E/eV 


Pz 


P<T 


Ptt 


d3 z 2 _ r 2 




dyz 


dxy 




0, 


Cu s 


Eci 


not. 


M 2 


1 


.47 





.053 

















.696 


.124 


.106 


.021 


Pa) 


S 


2* 


-2.40 





.273 





.039 











.323 


.083 


.119 


.163 


(Pad x 2_ v2 ) 


S 


3 


-2.94 








.659 











.341 














(P7rd X y) 


A 


4 


-2.96 


.704 











.245 

















.051 


(Pzd xz ) 


o 


5 


-3.76 





.086 


.003 


.344 











.057 


.006 


.017 


.487 


(d 3z 2_ r 2p CT ) 


S 


6 


-3.79 


.705 











.013 

















.282 


Pz 





7 


-4.11 





.635 


























.365 


Pa 


A 


8 


-4.32 








.933 


.001 














.049 





.017 


Ptt 


S 


9 


-5.93 

















.859 














.141 


dyz 


o 


10 


-6.17 








.199 











.801 














{dxyPir) 


A 


11 


-6.37 


.128 











.384 

















.488 


(d xz p z ) 


o 


12 


-6.44 





.390 





.409 











.104 


.035 


.062 





(P<jd 3z 2_ r 2) 


S 


13* 


-7.49 





.298 





.060 











.520 


.004 


.076 


.042 


(d I 2_j / 2Pff) 


S 
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TABLE 4: Comparison of experimental peak positions (in eV) with the LDA+U results at the high- 
symmetry points (ZRS and ZRT mean the Zhang- Rice singlet or triplet, respectively). 



r 



Orbital 


LDA+U 


Exp. 


{PnPa) 


-2.69 


-2.9 


Pz 


-3.83 


-3.9 


d x 2_ y 2 


-4.92 


-5.8 


dxy 


-5.40 


(PvPtt) 
(PvPtt) 


-5.57 


d( x ,y)z 


-5.87 


-6.5 


(7T,7r) 



Orbital 


LDA+U 


Exp. 


(P*d l x2 _ y2 ) (ZRS) 


0.65 


-1.2 


Ptt 


-2.43 


-2.4 


(Pzd {x . y)z ) 


-2.98 


-2.7 


{p-ir^xy ) 


-3.35 


(P*dl a _ ya ) (ZRT) 


-4.94 


-3.8 


(d( X .y) Z P Z ) 


-6.62 


-5.8 




-7.20 




-7.28 


-6.0 


(M) 


Orbital 


LDA+U 


Exp. 


(p*d x2 _ v2 ) (ZRS) 


-2.40 


-1.1 


iPzdxz) 


-2.96 


-2.5 


iP^dxy) 


-2.94 


-2.7 


Pz 


-3.79 


-3.8 


P-K 


-4.32 


Pa 


-4.11 


-3.8 


{dxyPir) 


-6.17 


-5.6 


dyz 


-5.93 


-6.6 


{dxzPz) 


-6.37 
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TABLE 5: TB-parameters obtained by fitting the LDA+U band structure and the VB photocmission 
spectra. The off-site energies in parentheses are the values from a fit only to the theoretical band structure 
in the cases where experimental corrections were approppriate. 



e 1 

t d 


b d 


£ d 


£D 




Eft 




£ P 


2.00 


-4.90 


-4.78 


-5.22 


-6.40 


-3.88 


-3.86 


-4.59 




>J 






1 7T7T 








1.33 


0.77 


0.71 


0.34 


0.37 
(0.32) 


0.84 
(0.77) 


1.15 
(0.77) 
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• Cu 

00 



LDA (both spin directions) 



2 - 



-2 - 



-4 - 



-6 - 



-8 - 



(0,0,1X0,0,0) 



(1,1,0) 



(1,0,0) 



(0,0,0) 



wave vector 



Pa 



Pa 




LDA + U (minority spin) 



2 - 




-6 - 



-8 - 



(0,0,1X0,0,0) (1,1,0) (1,0,0) (0,0,0) 

wave vector 



LDA + U (majority spin) 



- 




E? 

C 



-4 - 



-6 = 



-8 - 

(0,0,1X0,0,0) (1,1,0) (1,0,0) (0,0,0) 

wave vector 




Binding energy (eV) 



r / r - (7i,0) 

i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




1 " 1 1 1 1 1 1 1 " " " 1 1 1 1 1 r 

8 6 4 2 



Binding energy (eV) 



{%,%) I r -- {%,%) M 




1 (p dh ,) 

xi g x -y ' 

2 

3 (Pz d (x, y ,z) 
(P, d xy) 

4 (P d]y) 

5 ( d (x,y)zPz) 



6 (divP.) 



T 

8 



T 

6 



~r 

4 



~r 

2 



~r 





Binding energy(eV) 




Binding energy (eV) 




Binding energy (eV) 



'I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 




Binding energy (eV) 




Binding energy (eV) 




Binding energy (eV) 



Binding energy (eV) 




Binding energy (eV) 




